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We present recent results of SESAM's large scale lattice simulation of QCD with 
two dynamical flavours of Wilson fermions. 



1 Introduction 

As was pointed out by Gottlieb0, simulations with dynamical Wilson fermions 
at weak coupling on large lattices and at light quark masses are still lacking. 
In an effort to improve on this situation, SESAM (Sea-Quark Effects on Spec- 
trum and Matrix Elements) and T^L (Towards the Chiral Limit El) have been 
producing dynamical gauge configurations (with Wilson fermions) at a total 
of 4 sea-quark mass values, on two lattice volumes and working at a coupling 
which - as we shall show below - corresponds to a quenched coupling in the 
scaling regime. In this short note, w^-^iixnmarize some of SESAM's recent 



findings; more details can be found in 



2 Status of Simulation 

We work with two dynamical Wilson fermions on a lattice of dimensions 
16 3 x 32 and at a strong coupling of (3 — 5.6. SESAM is now close to com- 
pleting its QH2 (4x8x8 nodes) run-time of approx. 350 days and in this 
time we have produced, after thermalisation, 5000 trajectories of unit length 
(100 ± 20 molecular dynamics steps with dt — 0.01 - seetJ for details) at sea 
quark values jEsoa = {0.156,0.157,0.1575}. A half-time analysis was presented 
at Lattice 96B where we analysed {100, 160, 100} configurations per sea-quark 
taken at intervals of 25 units. This interval size is motivated from clean sig- 
nals in the autocorrelation function which emerge once the ensembles become 
larger than aprrox. 2000 trajectories. As an example we quote the integrated 
autocorrelation times of the plaquette rf n ^ = {3.3(5), 4.1(3), 7.1(5)} and - as 
a "worst case" - that of the average number of iterations in the HMC (using 
BiCGStab) if* = {20(1), 27(3), 31(4)}l. 
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3 Spectrum 

For each sea-quark mass we calculate hadron masses with valence quarks 
Kvai = {0.1555,0.1560,0.1565,0.1570,0.1575}. Gauge invariant smearing (50 
iterations and smearing parameter a=4) is applied to obtain smeared sources 
and smeared/local sinks. We use uncorrelated single-exp. fits to SL and 
SS data (simultaneously). Least-x 2 fits favour (linear + quadratic) chiral 
parametrizationscl for pseudoscalar and vector mesons as well as the nucleon 
(for vector particle and nucleon, a term oc m 3 / 2 does equally wellO). The lattice 





He am p a p 1 [GeV] Klight Kstrangc 


0.156 
0.157 
0.1575 


0.16065(8) 0.359(8) 2.14(5) 0.16058(8) 0.1576(2) 
0.15987(6) 0.341(8) 2.23(7) 0.15980(5) 0.1569(2) 
0.15963(11) 0.316(10) 2.44(8) 0.15944(7) 0.1544(1) 


quenched 




(3 = 6.0 


0.15718(6) 0.330(8) 2.33(6) 0.15709(4) 0.1544(1) 



Table 1: Lattice spacings and n- values. 



spacing is determined using the rho mass at n c ; using instead Kn g ht reduces 
a -1 by no more than 1%. The values for Ki; g ht and Estrange are obtained by 
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interpolating to the mass ratios —5- and — ^ and can be used to extract the 

quark massesE3. From table [j] we note that our lattice spacings correspond to a 
quenched f3 of around 6.0, one that is at the onset of the scaling regime. These 
values of a -1 are in agreement with those extracted from the interquark-forcdi 
An alternative way of taking the chiral limit is to use only data with 
K sca — K vai i n the extrapolation of mass ratios. We postpone this discussion 
until higher statistics are reached and an additional K sca from T^L becomes 
available. 

4 The 7r-Nucleon a term 

Over the past years there has been significant corunutational progress in the 
calculation of flavour singlet matrix elements (see til and references therein). 
It was shown, in particular, that the noisy-estimator technique with a random 
Zi noise is a promising method to calculate the amplitude of the disconnected 
contribution to oviV) where the term "disconnected" refers to : 

ovjv = m q (N\uu + dd-2ss\N) + 2m q (N\ss\N) 

connected , ^disconnected \ 

= ^TrN + Cr 7TJV ■ I 1 ) 
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Here, \N) is a nucleon state and m q denotes the current quark mass of the u 
and d (taken to be equal). 

Pioneering attempts to calculate^ the disconnected amplitude in full QCD 
have not led to overly clear signalsllj. Quenched Xhigh-statistics) calculations 
have produced rather more promising data (see c3 for a review and EjEJtJ) , 
however, it is unclear, a priori, what meaning can be attributed to disconnected 
amplitudes for gauge configurations which disregard the effects of dynamical 
quarks. In addition, quark masses are found to be much smaller in_full lattice 
QCD adding a significant uncertainty to the quenched calculations^]. 
Algorithmic Investigation 

We have carried out an extensive algorithmic investigation for the calculation of 
the disconnected diagram using our full QCD configurations to check whether 
- and if so, at what cost - a .decent signal can be obtained on a sample of 
maximally 200 configurations 0. The matrix element of the nucleon with qq 
insertions is obtained on the lattice from the following ratio : 

„, f ,_ 'L S {N(2,t)\Z z qq(z)\N(0)) _Ai„A-iM 
[) ~ £*Wz)|jV(0)> - 3rn' nA W 

^ constant + t (N\qq\N) , (2) 

where A(t) is the nucleon propagator. The hadron correlator is calculated in 
a standard fashion (using smearing as above). The disconnected contribution 
to the three-point function in the numerator is given by the correlation of the 
nucleon propagator and the disconnected insertion (Tr A xx ). In the follow- 
ing we present the results of our numerical investigation for the most efficient 
calculation of such disconnected diagrams. 

Noisy Estimator Techniques : The noisy estimator technique uses complex ran- 
dom sources with the property hirj.jv E ->oQ Vi(Ei ^)Vj(^> C) — to 
calculate (if M~ l rj) Ne times per configuration C. Using 157 configurations 
at our intermediate sea quark mass (i.e. ^f- = 0.76(1) and m q ~ 1.3 m s ) 
we have first compared Gaussian and Z2 noise sources to calculate the chiral 
condensate x ^ Tr(M -1 ). Monitoring the standard deviation of x versus the 
number of estimates Ne we found Z2 to outdo Gaussian noise by a near factor 
of 2; we pushed Ne all the way up to 300 in this investigation. Next, we varied 
the accuracy of the inversion residual r = which, obviously, we wish 

to relax as much as possible. Monitoring the quantity 6\ = xi r = 10~ 5 ) — x{ r ) 
we find that we can choose r ~ 10~ 4 and be safely within the 1 a error margin 
of the Z2 technique (Ne = 300). We now turn to the quality of the signal 
for lZ(t) which can be obtained with our 157 dynamical configurations and to 
the dependence of the signal-quality on the number of estimates applied; recall 
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Figure 1: R(t) measured with the noisy estimator technique (Ng = 300,) on 100 and 157 
gauge field configurations. Linear fits were performed in the range t = 7 to t = 15. The 
furthermost plot to the right shows a simulation with 157 quenched gauge-configurations 

((3 = 6.0 at K Ri Kstrangc/ 



that sa far, we have chosen Ne — 300, a number suggested by the authors 
of reflUj. In figure [I] we show lZ(t) measured on 100 and 157 configurations. 
With more than 100 configurations the signal becomes acceptable and the er- 
ror starts to display a behaviour. For comparison we show a plot from a 
quenched simulation with the same number of configurations (157); note that 
the signal is much worse, indicating that the determinant in the probability 
density with which configurations are sampled has a smoothing effect. Since 
a decent signal emerges, we monitor in figure || the chiral condensate and its 
variance as a function of Whereas the mean value of TZ(t) is unaffected 
when varying Ne , figure shows clearly that it does not pay to increase iV^ 
beyond 10-20 ! This drastically reduces the cost of the calculation. An analo- 
gous study using the volume source technique Ej, where M(C)x = (f> is solved 
for a volume source vector fa = 1, shows a much worse signal. Our method of 
choice is therefore the stochastic estimator technique with about 20 Zi sources 
per configuration and a relaxed residual r = 10~ 4 . The results of our simu- 
lation, where we calculate connected and disconnected amplitudes in full and 
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Figure 2: TrM 1 and iis variance as a function o/l/TVg (^57 gauge configurations). 



quenched QCD will be presented in a forthcoming paper. 
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